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ABSTRACT 



This paper shows that much new information about the dynamics of combat between 
two homogeneous forces modelled by Lanchester-type equations of modern warfare (also 
frequently referred to as "square-law 11 attrition equations) with temporal variations 
in fire effectivenesses (as expressed by the Lanchester attrition-rate coefficients) 
may be obtained by considering Liouville's normal form for the X and Y force-level 
equations. It is shown that the relative fire effectiveness of the two combatants 
and the intensity of combat are two key parameters determining the course of such 
Lanchester-type combat. New victory-prediction conditions that allow one to forecast 
the battle T s outcome without explicitly solving the deterministic combat equations 
and computing force-level trajectories are developed for fixed-force-ratio-breakpoint 
battles by considering Liouville's normal form. These general results are applied to 
two special cases of combat modelled with general power attrition-rate coefficients. 

A refinement of a previously known victory-prediction condition is given. Temporal 
variations in relative fire effectiveness play a central role in these victory-predic- 
tion results. Liouville's normal form is also shown to yield an approximation to 
the force-level trajectories in terms of elementary functions. 



1 . Introduction. 



Even though combat between two military forces is a complex random process ( see 
Note 1 of TAYLOR and BROWN^~^), as a consequence of the pioneering 1914 work of 

r ? f) l 

F. W. LANCHESTER 1 ( see Note 1) military operations analysts since about the end of 

World War II have used simplified deterministic differential equation models to develop 

insights into the dynamics of combat ( see Note 2). Today Lanchester-type models of 

quite complex military systems have been developed in the United States ( see , for example, 

[9] 

BONDER and HONIG ) and require a digital computer for their implementation. A simple 
combat model, however, may yield an understanding of important relations that are 
difficult to perceive in a more complex model, and such insights may provide guidance 

r o i 

for higher resolution computerized investigations ( see , for example, BONDER and FARREL 1 
and WEISS In this paper we will examine such an idealized Lanchester-type model 
in order to obtain some insights (specifically, the tradeoff between quality and quantity 
of weapon systems) into the dynamics of combat between two homogeneous forces with 
temporal variations in weapon system effectivenesses. 

In this paper we develop new victory-prediction conditions that sometimes allow 
us to forecast the battle’s outcome without explicitly solving the combat equations 
( see Note 3) and computing force-level trajectories ( see Note 4). We obtain these 
results for variable-coefficient Lanchester-type equations of modern warfare by consid- 
ering Liouville’s normal form of the X and Y force-level equations ( see p. 270 of 
INCE^~^ and KAMKE^^) and using techniques recently applied to Lanchester combat 
theory by TAYLOR and PARRY These results complement and extend those of Taylor 

and Parry , and they show that the key parameters affecting a battle’s outcome, at 
least for fixed-force-ratio-breakpoint battles with the initial force ratio held con- 
stant, are the relative effectiveness of the weapon systems and the intensity of combat. 

Such results are not only important in their own right but also useful in the quanti- 

[44 — 46 ] 

tative analysis of tactics ( see , for example, TAYLOR ). 



This paper is organized in the following fashion. First, we review F. W. 
Lanchester' s classic mathematical model of combat between two homogeneous forces and 
its extension to cases of time-varying fire effectivenesses. Next, we transform first 
the independent variable (time, t) and then the dependent variable (force level) to 
obtain Liouville ? s normal form for the X and Y force-level equations. Then we 
develop some new victory-prediction conditions from Liouville's normal form and apply 
these general results to two special cases of power attrition-rate coefficients. 
Finally, we discuss the significance of our developments. 

2 . Lanchester 's Classic Formulation . 

T 26 1 

F. W. Lanchester 1 hypothesized ( see Note 5) in 1914 ( see Note 6) that combat 
between two military forces "under modern conditions" could be modelled by ( see Note 7) 

dx/dt = -ay with x(t=0) = x^, 

( 1 ) 

dy/dt = -bx with y(t=0) = y^, 

where t = 0 denotes the time at which the battle begins, x(t) and y(t) denote 

the numbers of X and Y at time t, and a and b are nonnegative constants that 

are today called Lanchester attrition-rate coefficients and represent each side's fire 

effectiveness. Lanchester considered this simple model in order to provide insights 

into the dynamics of combat under "modern conditions" and justify the principle of 

concentration ( see Note 8). We will accordingly refer to (1) as Lanchester ? s equations 

of modern warfare . Various sets of physical circumstances have been hypothesized to 

yield them: for example, (A) both sides use aimed fire and target acquisition times 

[54] 

are constant ( see WEISS ), or (B) both sides use area fire and a constant density 
defense ( see BRACKNEY^^). Other forms of Lanchester-type equations have appeared in 
the literature, but we will not consider these here ( see DOLANSKY^^ and TAYLOR 
From (1) Lanchester deduced his famous square law 

b( x o “ x 2 (t ) } = a{y2-y 2 (t)}, (2) 



2 



which has the important implication that a side can significantly reduce its casualties 
by initially committing more forces to battle. It follows from (2) that 

x 0 

X will be annihilated — < 

y o 

Unfortunately, no simple relationship similar to (2) holds in general for variable 
attrition-rate coefficients so we consider other means for developing (3) . As is well 
known, the X force-level, x(t) , is given by 

x(t) = x n cosh(/ab t) - y n /^sinh(/atT t) . (4) 

U U > b 

We may also deduce (3) by writing (4) as 

x(t) = y{(x () -y 0v /|-)exp(v / ab t) + (x 0 +y o y^)exp(-/ab t)} , (5) 

and observing that x(t) can become zero if and only if the coefficient of the increas- 
ing exponential is negative. We observe from (5) that annihilation occurs in finite 
time. 

[53] 

As H. K. WEISS has emphasized, engagements that continue until one side is 

wiped out are rare. Thus, we see that a model of battle termination is required. 

Although we are well aware that battle termination is a complex random process for 

which it is by no means certain that force levels are the only significant variables 

( see Note 9), we assume that combat ends when either of two given "breakpoint" force 

ratios is reached. Introducing the force ratio u = x/y, we have that these "break- 

f f 

point" force ratios, denoted as u^. when X wins and u^ when Y wins, satisfy 
0 £ u^ < Uq = u(t=0) < u^ ^ +°o. Corresponding to a fight until the annihilation of 
one side or the other is the case in which u^ = 0 and u^ = -H 30 . 

Let us now consider such a fixed-force-ratio-breakpoint battle. Introducing the 
force ratio, u = x/y, we see that it satisfies the Riccati equation 

du/dt = bu 2 - a with u(t=0) = u Q = x^/y^. (6) 

Let u + = /a Jb denote the positive root of the quadratic equation bu 2 - a = 0 and 
observe that du/dt <0 <=> u < u + ( see Figure 2 of Taylor and Parry In 




3 



particular, => du/dt(t) ^ du/dt(t=0) < 0, whence follows 

THEOREM 1: Consider Lanchester's equations of modern warfare with constant 

attrition-rate coefficients (1) . Then X will lose a fixed-force-ratio 
breakpoint battle in finite time if and only if x q/Yq k ^a/b . 

We observe that (3) is a special case of Theorem 1 corresponding to u^ = 0 and 
u^ = +°°. One result of this paper is to generalize [in a way different from that given 
by Taylor and Parry ( see Section 6)] Theorem 1 to cases of variable attrition-rate 

coefficients . 

3. Variable Attrition-Rate Coefficients . 

[3 5 81 

The pioneering work of S. BONDER 1 ’ ’ on methodology for the evaluation of 

military systems (especially mobile systems such as tanks, mechanized infantry combat 
vehicles, etc.) has generated interest in variable-coefficient Lanchester-type equations 
and has led to improved operations research techniques for the prediction of these 
coefficients ( see Note 10). Let us therefore consider 

dx/dt = -a(t)y with x(t=0) = x , 

U (7) 

dy/dt = -b(t)x with y(t=0) = y^, 

where a(t) and b(t) denote time-dependent Lanchester attrition-rate coefficients. 

These coefficients depend on such variables as force separation, tactical posture of 
targets, rate of targer acquisition, firing doctrine, firing rate, etc. ( see reference 8). 
We will also refer to (7) as the equations for a square-law attrition process , since 
an "instantaneous" square law holds even when a(t)/b(t) is not constant ( see Taylor 
and Parry ; also references 45 and 47). 

A large class of tactical situations of interest can be modelled with the follow- 
ing general power attrition-rate coefficients ( see reference 8) 

a(t) = k (t+C) P , and b(t) = k, (t+C+A) V , (8) 

a b 
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where A,C ^ 0. We will call A the offset parameter , since it allows us to model 
(with p,v ^ 0 ) battles between weapon systems with different effective ranges ( see 
Note 11). We will call C the starting parameter , since it allows us to model (again, 
with p,v ^ 0 ) battles that begin within the maximum effective ranges of the two systems. 
Restrictions that must be placed on p and v, which are not necessarily integers, 
are discussed below. 

The above nomenclature is motivated and possible applications of our results are 

[3 5 ] 

indicated by considering S. Bonder’s * model of the constant-speed attack on a 
static defensive position 

dx/dt = -o(r)y = -a 0 (l-r/R a ) P y, dy/dt = -B(r)x = -B^l-r/R^x, (9) 

where r denotes the range between opposing forces, P,v ^ 0, and denotes the 

maximum effective range of the Y weapon system. Range is related to time by 
r(t) = Rq - vt , where R^ denotes the opening range of battle and v > 0 denotes the 
constant attack speed. Then the offset and starting parameters are given by 

A = (Rg-R a )/v, and C = (R^-R^/v. (10) 

We observe that A,C £ 0 if and only if R 0 ^ R £ R~. By considering (10) and 

p cc u 

Figure 1, the reader should have no trouble in understanding our terminology for A 
and C. In the model (9) p, for example, is used to model the range dependence of 
Y’s attrition-rate coefficient ( see Figure 2). 

From (7) we obtain the X force-level equation 

^7 ~ Ojj~£n a(t)} , ~- a(t)b(t)x = 0 , ( 11 ) 

with initial conditions 

x(t= 0 ) = x Q and { [ 1 /a ( t) ]dx/dt } t= Q ~ “Yq» 

where t^ = max(t^,tQ) , and t^ denotes the right most finite singularity of the 
X force-level equation ( see Ince^^). The coefficients a(t) and b(t) are then 
positive continuous functions V t > t^. For the coefficient ( 8 ), we observe that 



5 
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Figure 1. Explanation of offset parameter A and starting parameter C for power attrition-rate coefficients 
modelling constant-speed attack. [Notes: 1. The maximum effective ranges of the two weapon 
systems are denoted as R^ and R 3 . 2. The opening range of battle is denoted as Rq and 
(as shown) Rq < minimum (R a ,Rg). 3. The offset parameter is given by A = (Rg-R^/v. 4. The 
starting parameter is given by C = (R a -RQ)/v.] 
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imum effective range of the system is denoted as R a = 2000 meters. 2. a(r=0) = ag = 0.6 X 

ualties/ (unit time x number of Y units) denotes the Y force weapon system kill rate at 
o force separation (range). 3. The opening range of battle is denoted as Rg = 1250 meters 
(as shown) Rq < R a .] 



tg = -C. Moreover, to insure the existence of all derivatives required in subsequent 
analysis, we assume that the second derivatives of a(t) and b(t) exist V t > t^. 
We similarly have the Y force-level equation 

g-{^lnb(t)}^-a(t)b(t)y-0, (12) 

with initial conditions 



y (t=0) = y Q and { [1/b (t) ] dy/dt > t=0 = -x Q . 



It is necessary to place further restrictions on a(t) and b(t) in order to 
insure that, for example, the transformation introduced in the next section is well 
defined. Thus, we assume that the following condition holds. 



CONDITION (A): 



t 

* 

a(s)ds and 



t 

b(s)ds are bounded for all finite t £ t^. 



[221 

This condition also guarantees ( see Theorem 6.4.2 on p. 226 of HILLE ) that (7) has 
a continuous solution for all finite t £ 0 ^ t^ ( see Note 12) . If Condition (A) is 
to hold, then for the general power attrition-rate coefficients we must have y,v > -1. 



4. Transformation of the Battle’s Time Scale. 



Let 



t = 



/a(s)b(s) ds , . 



(13) 



and denote T(t=0) as Tq. It follows that Tq ^ 0 for t^ £ 0. Condition (A) 

guarantees that t = x(t) is well defined (i.e. bounded for all finite t ^ t^) by 

[21 

the Cauchy-Schwarz inequality for integrals ( see p. 123 of BELLMAN ). The transfor- 
mation is invertible, since dx/dt = /a (t )b (t ) > 0 for t > t^. Thus, we may consider 
that t = t(x). We assume (and give below conditions that guarantee) that 
lim x(t) = -foo so that the range of x(t) is [0,-H») for te [t n ,-H») . 

Considering the constant coefficient result (4) , we will call the quantity 
A(t)b(t) the " intensity of combat" ( see also Taylor and Parry^^) ; since the larger 
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it is, the more quickly the battle is moving towards termination. Then T ” t q is 

related to the average intensity of combat by 

t 



T - T Q = { (1/t) 



/a(s)b(s) ds } t = v / a(t)b(t) t. 



(14) 



Applying the transformation (13) to (11) and (12) , we obtain 



d z x rid. 

-7— 7 + 1 ~r~ &n 

dx z 2 dx 



)?-* - 0 . 

_a(t)J dx 



(15) 



with initial conditions 



x(x=x Q ) = x Q and { [b(t)/a(t) ] 1/2 dx/dx} 



T = T 



= -y 



o* 



and 



A , f l d ra(t)Td Z . „ 
d^ +{ 2 d7* n = °’ 



with initial conditions 



, 1/2 



(16) 



(17) 



y(T=x Q ) = y Q and {[a(t)/b(t)] dy/dx}^^ = -x Q . 

Taylor and Brown h ave shown that the X force-level equation (11) may be 
transformed into a linear second order differential equation with constant coeffi- 
cients if and only if 

- r-tn - CONSTANT, 

/aYt)b(C) dt L b<C) J 

with the desired transformation being given by t = K j /a (s)b (s) ds , where 
denotes an indefinite integral and K an arbitrary constant. Hence, equation (11) 
may be transformed into such a constant coefficient equation if and only if (15) is a 
constant coefficient equation. We observe that (15) has a solution in terms of ele- 
mentary transcendental functions in the quasi-autonomous case in which the ratio of 
attrition-rate coefficients is constant (see Note 13) , i.e. 



ft 



ds 



b(t)/a(t) = CONSTANT. 



(18) 
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Except when (17) holds, the solution to (11) is complex, and the qualitative behavior 
of force-level trajectories has been difficult to establish. In particular, one is 
interested in answering such questions as 

(Ql) Who will "win"? Be annihilated? 

(Q2) How do force levels decrease over time and how many survivors will the 
winner have? 

(Q3) How do changes in the initial force levels and/or weapon system parameters 
affect the outcome? Is concentration of forces a good tactic? 

(Q4) How long will the battle last? 

We will now show how to answer question (Ql) without explicitly solving the Lanchester- 
type equations (7). 

Equation (15) is highly significant because it clearly shows that the course 
of combat depends on just the two weapon system parameters: (1) R(t) = a(t)/b(t), 
the relative effectiveness (Y to X) of the two weapon systems, and (2) I(t) - 
/a(t)b (t) , the intensity of combat (through equation (13), which relates I(t) to 
x) . Both these parameters may vary over time, and equation (15) tells us that the 
nature of such temporal variations in relative effectiveness has a significant effect 
upon the course of combat. Moreover, this relationship may be more explicitly seen 
by transforming (15) to Liouville’s normal form. 

5 . Reduction to Liouville T s Normal Form . 

Let us assume that Tq > 0 (i.e. 0 > t^) so that a(t)/b(t) is twice differ- 
entiable and satisfies 0 < a(t)/b(t) < +°° for 0 £ t < +°°. Let a^ denote a(t=0), 
etc., and recall that t = t(x). The substitution 




(19) 



r 

transforms (15) into the so-called normal form ( see RAINVILLE l j ) with the first 



derivative of the dependent variable removed 




( 20 ) 
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with initial conditions 



X(T=T 0 ) = X 0’ and dX/dx(x=x Q ) = -y 0 /a 0 /b 0 -x 0 e 0 . 



where 



F(t) = P"(t)/P(t), 
R(t) = a(t)/b(t) , 



-1/4 



e(t) = 



P(t) = [R(t) ] 

1 d 



4/a(t)b(t) 



dt 



2.n R , 



( 21 ) 

( 22 ) 



e Q denotes e(t-O), and P'(t) denotes dP/dx. Equation (20) is Liouville's 



celebrated normal form ( see p. 23 of LIOUVILLE 1 J ). Similarly, 

-1 1/4 

D(,U;/D 

y(x) = y ( x) 



a(t)/c 



0 



yields 



with initial conditions 



d 2 Y 



dx^ 



- {1 + G(x)}Y = 0, 



(23) 



(24) 



where 



Y ( t=t 0> = y 0’ and dY/dx (x=x Q ) = -x 0 ,/b 0 /a 0 + y 0 e 0’ 



G(x) = Q"(x)/Q(x), and Q(x) = [R(t)] 1/4 = 1/P(x). 



(25) 



We observe that 



F(t) + G(x) = 2(7-lnP) 2 . 

dx 



(26) 



It is sometimes convenient to express F and G in terms of the old time 
variable t. Then 



and 






G(x) = o_V zkm \ { jr^n R (t ) } 2 - F (x) . 



8a(t)b(t) dt 

Observing that (21) may also be written as 



F(x) = -j-? £n P(x) +{ -j- £n P(x) } 2 , 
dx dx 



(27) 



(28) 



(29) 
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we see that d £n P/dT = constant (or, equivalently, e(t) = constant) implies that 
F(t) is constant, although the converse may not be true. 

Writing (20) as d X/dT -X = F(t)X, we may use variation of parameters ( see 

r 23 1 

pp. 122-123 of Ince L J ) to obtain the solution to (20) as 



X(x) = X Q cosh (t-Tq) - [y 0 t/a 0 / b 0 + x o £ 0 l sinh(T_ 'V + 



F(a)sinh(T-a)x(a)da. (30) 



0 

In terms of the original time variable t and dependent variable x, we have 
-1 1/4 



x(t) = 



a(t) /a 



b(t)/b 



{xq cosh ( 



[{ 77 ^ 7 - [ V )sinh( /a(o)b(a) da) x(s)ds. (31) 

[b (t)J J ds a(s)b(s) ds |_a(s)J J 



t /a ft 

/a(s)b(s) ds) - [y Q /^-+ x 0 e Q ]sinh( /a(s)b(s) ds) } 

0 ° 0 

1 1/4 



The value of the Volterra integral equation (30), however, is not so much for 

direct computation as it is for suggesting an approximation, the so-called Liouville- 

[321 

Green approximation ( see OLVER ), to the solution of the X force-level equation 

(11). If the appropriate fractional power of the relative effectiveness is "slowly 

varying," then by (21) we would expect that |f(t)| << 1 so that we could drop the 
integral terra in (30) to obtain 

X(T) = x 0 cosh (T-T 0 ) - [y 0 ^7b^ + x 0 e 0 ]sinh(T-T 0 ), (32) 

where X(t) denotes the Liouville-Green approximation. A theoretical error analysis 
[321 r 33 1 

appears in Olver ( see also OLVER ), although we will not pursue this matter 

further here. We observe that F(x) ^ 0 V x ^ Xq implies that while X(x) ^ 0 we 

have X(t) ^ X(t) and similarly when F(t) £ 0. As we shall see below, such cases 
in which F(x) is always ^0 or £ 0 are readily encountered in applications. In 
terms of the original variables t and x, we have 

+ x 0 e 0 ]sinh(j /a(s)b(s) ds)} , (33) 

0 



x(t) = 



a(t) /a 
b(t)/b 



{Xq cosh ( 



/a(s)b (s) ds) “ [y n / 

0 
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which we may also write in terms of the average combat intensity, /a ( t)b(t ) = 
ft 



(1/t) 
x(t) = 



/a (s)b (s) ds , as 

0 r / w nl 
a (t ) /a 



0 



b(t)/b 



0 



{xq cosh (/a(t)b (t) t) - [ y Q / x Q e 0 ]sinh (^ a (t) b (t) t)} . 



( 34 ) 



6. Some New Victory-Prediction Conditions . 

In this section we show that much valuable information about the course of 
combat (for example, force-annihilation prediction) may be obtained directly from 
Liouville's normal form (20) without making any kind of approximation. Let us assume 
that the model (7) holds for all time and consider a fixed-force-ratio-breakpoint 
battle ( see Section 2). In order to insure that the battle terminates in finite time, 
we must make certain technical assumptions about the mathematical nature of the attri- 
tion-rate coefficients in (7). For the reader's convenience, we list in Table I the 
principal such conditions that we use in subsequent developments. We observe that all 
these conditions hold for constant attrition-rate coefficients. Let us also observe 
that Conditions (B) and (ND) imply Condition (C) . 

[51 1 

We record here for future reference Taylor and Parry’s J generalization of 
Theorem 1 to cases of variable attrition-rate coefficients ( see Note 14). 



THEOREM 2: Assume that Conditions (B) and (ND) hold. Then x^/y^ < ^ a o^O 

implies that the X force will lose a fixed-force-ratio-breakpoint battle 
in finite time. 



PROOF: Introducing the force ratio, u = x/y, we have 

du/dt = b(t)u 2 -a(t) with u(t=0) = u^ = x^/y^. (33) 

Let u + (t) = /a (t ) /b (t) = /R(t) denote the positive root of the quadratic equation 
b(t)u 2 -a(t) = 0, and observe that du/dt < 0 for any u < u + (t) ( see Figure 2 of 
Taylor and Parry ^^). Condition (ND) means that u_^(t) is nondecreasing. It is 
readily shown that du/dt (t r -0) < 0 and u_^_(t) nondecreasing imply that du/dt (t) < 0 
for all t £ 0 ( see pp. 526-527 of Taylor and Parry ^ 51 3- Consequently, from (35) we 
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TABLE I. 



Assumptions about Attrition-Rate Coefficients. 



CONDITION (A) : 



CONDITION (B): 



CONDITION (C) : 



CONDITION (ND) : 



t t 

a(s)ds and j b(s)ds 



for all finite t ^ t^. 



lim 

t->+°° 



b(s)ds = +». 



lim 

t-H<» 



/a(s)b(s) ds = -H». 



R(t) is nondecreasing on 



are bounded 



[ 0 ,+~) . 
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see that x Q /y 0 < /a Q /b 0 and Condition (ND) yield that du/dt(t) < 0 V t ^ 0. It 
then remains to be shown that X's breakpoint force ratio is reached in finite time 
( see Note 15). Let us observe that R(t) nondecreasing and 0 < a(t),b(t) < +°° V t (0,+°°) 
» a^ < +“ and b Q > 0. From Condition (ND) and the fact that du/dt(t) < 0V t 0, 
we have 



^r(t) = b(t){u 2 -R(t)} 






0“0 a o } 



b(t) du 



0 



dt 



(t=0) 



Hence 



so that 



u(t) " u 0 + 1 (df) dt s U 0 + ^r3? (t ' 0> 



ft 



lim 

t-H-® 



b (s)ds = + 00 implies that u(t) goes to 



r t 

b (s)ds , 

0 

^ 0 in finite time. 

Q.E.D . 



We observe that Theorem 2 says that the X force will be annihilated in finite 

XX f 

time (i.e. u(t=t ) = 0 for t finite when u = 0) . 

na Y 

We will now deduce some new results that many times allow us to predict such a 
battle’s outcome from the initial conditions. Let U = -(l/X)dX/dr so that (20) 
becomes the Riccati equation 

+ e Q . (36) 

We observe that 



^ = U 2 - 1 - F (x) 
dx 



with U(t=t ) = — / 

0 V b o 



U(x) = ^ R 1/2 + e(x) . (37) 

If F(x) £ 0 V x 2: x , then dU/dx £: U 2 -l, and U(x=Xg) = Uq > 1 » lim U(x) = +“ 

T->T 

X 

a 

for finite. We observe from (22) and the definition of t^ that R(t(T)) and 

a 

e(t) are finite for all t ^ > 0 so that by (37) the X force must be annihilated 

in finite time, since by Conditions (A) and (C) T(t) is a strictly increasing func- 
tion V t e [tg,+°°) with range [0,+°°). If we additionally assume that dR/dt ^ 0, 
then du/dt (t ) < 0 for all t e[0,t^] with lim u(t) = 0 so that the X force will 

a 
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lose any fixed-force-ratio-breakpoint battle in finite time. The proof is by 
contradiction: dR/dt ^ 0 => u + (t) nonincreasing => if there exists t q such that 
dx/dt(t^) ^ 0, then du/dt(t) ^ 0 V t ^ t^ impossible to have u -> 0. Thus, we 
have proved 



THEOREM 3: Assume that Conditions (A) and (C) hold. If F(x) ^ 0 V t ^ x , 

x o F 

then — (1-e ) < / — implies that the X force will be annihilated in finite 

y o 0 / b 0 

time. Furthermore, if dR/dt ^ 0 V T ^ x , then the X force will lose any 



fixed-force-ratio-breakpoint battle in finite time. 



If F(t) £ 0 and dR/dt i OVt i T n , then dU/dr £ U 2 - 1 and U„ < 1 => 3 t 

U u 

finite such that U(t^) = 0 with U(x) > 0 V x e [ T Q> T y)- Since dR/dt ^ 0 =* e(t) ^ 0, 
we see from (37) that Y must be annihilated in finite time, since x(t) is a one- 
to-one mapping of [t^,+°°) onto [0,-H»). By observing that 



du , \ 9 r 1 dR 

dF = b(t)u { 2l(o^dr 



+ 1-U 2 + (-p-£n P) 2 }, 
dx 



we see that du/dt(t) > 0 V t ^ 0 as long as u > 0. Hence u(t) is strictly increasing 

3i 3 . 

for te [0,t Y > and lim u(t) = -4^° for t finite so that Y will lose any battle 

c-*t- 

f a 

with u > u n . Thus, we have proved 
X U 



THEOREM 4: Assume that Conditions (A), (B) , and (ND) hold. If F(x) ^ 0 V x^x^, 

x /a 

then — (1-e^) > /: — implies that the Y force will lose a fixed-force-ratio- 

y o 0 >' b o 

breakpoint battle in finite time. 



By considering Lioville’s normal form (24) for the Y force-level equation, 
one may similarly prove Theorems 5 and 6. 



THEOREM 5: Assume that Conditions (A) and (C) hold. If G(x) £ 0 V x^ x^, 

then — > (1+e^) implies that the Y force will be annihilated in finite 

y o v b 0 0 

time. Furthermore, if dR/dt ^ 0 V x^ x^ (i.e. Condition (ND) holds), then Y 



will lose any fixed-force-ratio-breakpoint battle in finite time. 
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a3 >h 



t 



THEOREM 6: Assume that lim 

f-H-co 



a(s)ds - +°° and that Condition (A) holds. 



0 



x / a 

G(t) ^ 0 and dR/dt ^ 0 V T ^ t q , then — < y — (l+e Q ) implies that the 
force will lose a fixed-force-ratio-breakpoint battle in finite time. 



If 

X 



Although F(t) and G(t) are related by (29), for all the attrition-rate 
coefficients that we have so far considered, F(t) £ 0 G(t) £ 0. Also, in all these 
cases F(t) ^ 0 « dR/dt ^ 0. Since x Q /y 0 < /a Q /b 0 (l+c Q ) =» (x Q /y Q ) (l-c Q ) < /a Q /b 0 
and (x^/y^) (1-e^) > /a^/b^ » x^/y 0 > /a^/b^ (1 +Eq) , we see that Theorem 6 is stronger 
than Theorem 3 and Theorem 4 stronger than Theorem 5. Under these conditions, Theorems 
4 and 6 provide no additional information on battle outcome over that contained in 
Theorems 3 and 5. 



7. Application to General Power Attrition-Rate Coefficients . 

We now apply the above results to combat between two homogeneous forces modelled 
by (7) with the general power attrition-rate coefficients (8) . We distinguish between 
two cases: (I) power attrition-rate coefficients with no offset (i.e. A = 0) , and 

(II) power attrition-rate coefficients with the same parity (i.e. p = v) and positive 
offset. In order to invoke Theorems 3 through 6, we must have C > 0 for the general 
power attrition-rate coefficients (8) (cf. (19)). In order that Condition (A) be 
satisfied, we must have p,v > -1; and then Conditions (B) and (C) hold ( see Table I). 

7.1. Power Attrition-Rate Coefficients with No Offse t . 

In this case we have 

a(t) = k (t+C)^ and b(t) = k^(t-HC) V , (38) 

so that 

dR/dt ^ 0 p ^ v. (39) 



Theorem 2 then yields 



COROLLARY 2.1: For the power attrition-rate coefficients (38) with p ^ v and 

x 

C > 0, — < /-^ implies chat Che Y force will win a fixed-force-ratio- 

*0 Vb o 
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breakpoint battle in finite time. In particular, the X force will be 

annihilated in finite time when = 0. 

In preparation for invoking Theorems 3 through 6, we compute 

_ (y-v) (3y+v+4) 

F(T) ■ 4(y+v+2 )*T* ’ 



and 



r , T s _ (v-u) (y+3v+4) 
G(T) ” 4(y+v+2)^T‘ i * 

since (13) with (38) yields 



T - T(t) -4!^ (t + C) <U+V+2)/2 , 



(p+V+2) 



,(y+v+2)/2 



with = (2/k^k^/ (y+v+2) )C “ > 0. We observe that 

F (t) £ 0 and G(t) ^ OVt > 0 « 



(40) 



(41) 



(42) 



(43) 



Thus, as discussed above at the end of Section 6, Theorems 4 and 6 provide no additional 
information on force annihilation over that contained in Theorems 3 and 5. Furthermore, 
Theorem 3 is stronger than Theorem 2, since shich is given by (22) evaluated at 

t = 0, is ^ 0. Hence, we omit any corollary to Theorem 3. As a corollary to 
Theorem 5 we have 



COROLLARY 5.1: For the power attrition-rate coefficients (38) with y £ v 

and C 



-L- ^ ♦ J- • i- Wi. UUU UL LL J.LXV11 I- U. LL, i. XU LO \ W / W J. L I 1 \JL • 

> °’ ~ > t c~ (M+V+2)/2 } implies that the X fore 



y 0 ' “0 

will win a fixed-force-ratio-breakpoint battle in finite time. In particular. 



4» / k^k^ 



the Y force will be annihilated in finite time when u = +°°. 

A 



Let us also write the X-victory-prediction condition of Corollary 5.1 as 

f(C). (44) 



X 0 . / a c (y-v)/2 ( (y-v) c ~(v+l) = 

y Q V kfe 4 k, 



4k b 



Observing that lim f(C) = +°° for y > v, we see that the victory-prediction condi- 
C+0+ 

tions of Corollaries 2.1 and 5.1 become stronger as C decreases and are meaningless 
for C = 0. 
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We now show that other considerations, however, readily yield force-annihilation 
conditions when C = 0. [Unfortunately, without further analysis these results are 
limited to force— annihilation prediction and do not apply directly to fixed— force— ratio- 
breakpoint battles with u^ < +°° and u^ > 0.] The solution to (7) with the power 
attrition-rate coefficients (38) may be written for C = 0 as ( see Note 16) 

x(t) = (p P /(2p)){C x -l-C Y }A g (T) + (p P /(2v/^)){C x -C Y }B B (T) , (45) 

where and denote the generalized Airy functions of the first and second kinds 

of order 6 ( see SWANSON and HEADLEY^ 41 !), p = (y+1) / (y+v+2 ) , 6 = (v-y)/(y+l), 

C x = x 0 r(l-p) , C Y = r(p)(A a k b /(y+xH-2)) 1_2p , and T = ( /Tk^/ (y+1) ) 2p t y+1 . 

r 41 1 

Observing that A y (5) , B (?) > 0V K * 0, lim A ( O = 0, and lim B (?) = +~, 

V £->+o° v 

we see from (45) that lim x(t) = -°° if and only if the coefficient of B (T) is 

t^+°° 

negative (i.e. C < C y ) . Hence, we may conclude (cf. the development of (3) from (5)) 



PROPOSITION 1: For a f ight-to-the-f inish modelled with the power attrition- 



rate coefficients (38) and C = 0, the X force will be annihilated in finite 

J,(£L 




[50] 



r(i-p) * 

For C > 0, (45) does not take nearly such a convenient form ( see Taylor and Brown lJU] ) , 

and the analogue of Proposition 1 for C > 0 would involve the generalized Airy func- 
tions and their derivatives evaluated at a positive argument. Thus, the prediction of 
force annihilation by this approach for C > 0 would require tabulations of higher 
transcendental functions. Such tabulations do not currently exist ( see reference 50). 
Thus, we see the usefulness of the "strong" annihilation conditions given by Corollaries 
2.1 and 5.1, which contain only elementary functions: "exact" annihilation-prediction 

conditions involve higher transcendental functions. 



7.2. Power Attrition-Rate Coefficients with the Same Parity and Positive Offset . 
In this case we have 

a(t) = k a (t+C)^ and b(t) = k^ ( t+C+A) ^ , (46) 



with A,C > 0. It follows that 
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dR/dt £ 0 « p £ 0. 



(47) 



Then Theorem 2 yields 



COROLLARY 2. 
p ^ 0 and 



2: For the offset power attrition-rate coefficients (46) with 

x /a 

A,C > 0, — < Jtt- implies that the Y force will win a fixed- 

Y r y 



r 0 



0 



force-ratio-breakpoint battle in finite time. 



Let us now see what information is yielded by Theorems 3 through 6. It is more 
convenient to express F and G in terms of the old time variable t ( see (27) and 
(28)). Then 

v = yA{4(y+2)(t+C) + (3u+4)A} 

u +2 u +2 ’ 

16k a k b ( t+C ) M ( t+C+A ) p A 

and 



g(t) = -pA{4(y+2) (t+C) + (y+4)A} _ 
1 6^1^ ( t+C ) M+2 ( t+C+A) P+2 

We observe that 



(49) 



F(t) ^ 0 and G(t) ^ OVt > 0 * p £ 0, (50) 

so that Theorems 4 and 6 again provide no additional information over that contained 
in Theorems 3 and 5. Furthermore, we may omit consideration of Theorem 3, since it is 
implied by Theorem 2. As a corollary to Theorem 5 we have 

COROLLARY 5.2: For the offset power attrition-rate coefficients (46) with 

x !<3l 

p ^ 0 and A,C > 0, ~ > \Z\T~ {H — } implies that the X 

y 0 0 4/a b " C (C+A) 

force will win a fixed-force-ratio-breakpoint battle in finite time. 

Let us also write the X-victory-prediction condition of Corollary 5.2 as 

x ys (14 - a/C )-^- ^ - ~ t . gw. <5i) 

We observe that lim g(C) = +°°, and we again see that the victory -prediction conditions 
C-*0 + 

of Corollaries 2.2 and 5.2 become stronger as C decreases and become meaningless for 
C = 0. 
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8. Discussion. 



ron 

In his classic 1914 paper , Lanchester assumed that the combatants’ fire 

effectivenesses (as expressed by the Lanchester attrition-rate coefficients) were 
constant over time and deduced his famous square law (2), which allows one to tradeoff 
quality versus quantity of weapon systems by means of the condition for equality of 
’’fighting strengths” 

x o /y o = (52) 

where a and b denote constant attrition-rate coefficients ( see Section 2). Thus, 
we see that equality of Lanchester-type fighting strengths ( see Note 17) depends on 
two parameters: the initial force ratio and the relative effectiveness. When the 

timing of military actions is considered, we add a third parameter, the intensity of 
combat, to this list of significant combat parameters. In the paper at hand, we 
extended these well-known constant-coefficient results to battles between two homogeneous 
forces with temporal variations in the fire effectivenesses. 

No such simple relationship like the square law (2) , which yielded (52) , holds 
in general for variable attrition-rate coefficients ( see Note 18). By transforming 
the independent variable to normalize the battle’s time scale by the intensity of com- 
bat, we found that the course of combat depends on two weapon system parameters: (I) 
relative fire effectiveness, R(t) = a(t)/b(t), and (II) intensity of combat, 1 ( t ) = 

/a^( t)b (t) . Moreover, when the temporal variations in relative fire effectiveness 
R(t) follow a regular pattern (e.g. R(t) nondecreasing), the battle’s outcome can 
many times be predicted from the battle’s initial conditions. To obtain such results, 
we considered Liouville’s normal form for the, for example, X force-level equation 
and found that it not only yields new battle-outcome-prediction conditions but also 
suggests an approximation to the time history of the X force level x(t) . As seen 
from ( 21 ), Liouville’s normal form ( 20 ) introduces second order conditions (e.g. the 
second derivative of relative fire effectiveness with respect to transformed time) into 
Lanchester-type combat analysis. 
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The new victory-prediction results ( see Theorems 3 through 6) that we have 
developed here are complementary to those of Taylor and Parry ( see Theorem 2). This 

complementary nature is seen by observing that under the appropriate comparable condi- 
tions we have 



from Taylor and Parry’s Theorem 2 (— £ 0): Y will win if 

at 

from our new results. Theorem 5 (G(x) £0): X will win if 




where 




(53) 



For both special cases of general power attrition-rate coefficients considered in 
Section 7 above we had dR/dt ^ 0 « F(x) £ 0 ** G(x) £ 0 ** ^ 0. Although these if- 

and-only-if statements do not hold in general, they do hold for these particular coeffi- 
cients. In both cases, we observe that for 






(54) 



we cannot say by this approach who will be the loser of the fixed-force-ratio-break- 
point battle ( see Figure 3) . We observe that force annihilation is a special case of 
these outcome-prediction conditions. From both (54) and Figure 3, we see that there 
is a "gap" in these victory-prediction conditions (i.e. Theorems 2 through 6). The 
price of removing this "gap," however, is the introduction of higher transcendental 
functions ( see Section 7 above and Taylor and Brown ) . Furthermore, "exact" results 
with no such gap in the victory-prediction conditions are apparently only possible 
for a f ight-to-the-f inish in which one side or the other is to be annihilated. 

We also refined a victory-prediction result of Taylor and Parry (as applied 
to the model (7) under study) by adding a restriction on the attrition-rate-coefficients, 
i.e. our Condition (B) , to the assumptions of Theorem 2. Taylor and Parry did not 
note that such a condition must be assumed in order to insure that the battle will 
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Figure 3. "Gap" in the Victory-! 



terminate. All our new results, i.e. Theorems 3 through 6, contain some such restriction 



on the attrition-rate coefficients. 

The results of this paper may be used in parametric analyses ( see BONDER 1 ' J ) of 
the dynamic combat interactions between two homogeneous forces with time-(or range-) 



,17], 



dependent weapon system capabilities. Such models are of particular interest in light 

of the work of S. BONDER^’^ and others ^ ^ 5 ^ > 37 ,38] the prediction of Lanchester 

attrition-rate coefficients from weapon system performance data and the work of 

[ 12 ] 

G. CLARK on the estimation of such (time-dependent) coefficients from Monte Carlo 

simulation output. A further discussion of applications is to found in references 8 
and 43. As is always the case, however, the insights gained into combat dynamics from 
such Lanchester-type models are no more valid than the models themselves. 



9 . Summary . 

In this paper we have developed insights into the dynamics of Lanchester-type 
combat between two homogeneous forces with temporal variations in weapon system 
effectivenesses by considering Liouville’s normal form for the X and Y force-level 
equations. Our principal results were some new outcome-prediction results that com- 
plemented those of Taylor and Parry We also saw that Taylor and Parry’s victory- 

prediction result (as applied to our model (7)), Theorem 2, must be refined by making 
an additional assumption about the attrition-rate coefficients. Liouville’s normal 
form also suggests an approximation, the Liouville-Green-Lanchester approximation, to 
the force-level trajectories. Additionally, by transforming the battle’s time scale, 
we saw that the relative fire effectiveness of the two combatants and the intensity 
of combat were two key parameters affecting the course of battle. Our new victory- 
prediction conditions for fixed-force-ratio-breakpoint battles were applied to two 
special cases of combat modelled with general power attrition-rate coefficients: (I) 

power attrition-rate coefficients with no "offset" (modelling, for example, two weapon 
systems with the same miximum effective range), and (II) offset power attrition-rate 
coefficients with the same parity (modelling, for example, weapon systems with different 
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